The Narrow Line Region of the Seyfert 2 galaxy Mrk 78. 

An infrared view 
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ABSTRACT 

We report near- infrared spectroscopic data for the Seyfert 2 galaxy Mrk 78, taken with the 
LIRIS near- infrared camera/spectrometer at the William Herschel Telescope (WHT). The long- 
slit spectra clearly show extended emission. The resolution and depth of the near- infrared spectra 
allows the examination of its morphology and ionization regions, and a direct comparison with 
similarly deep visible spectra. The emission-line ratios obtained are used to derive the extinction 
towards the nucleus. The detection of strong features such as [Fe II], H2, hydrogen recombination 
lines and the coronal [Si VI]A1.962 line is used to study the kinematics and excitation mechanisms 
occurring in Mrk 78, revealing that despite of the strong radio-jet interaction present in this 
object, photoionization from the active nucleus dominates the narrow line region emission, while 
UV fluorescence is the source of the H2 emission. Lines with extended emission yield velocity 
distributions with an amplitude of about 600 km s _1 , the consequence of an eastern lobe moving 
away from us plus a western lobe with the opposite contribution. We used the photoionization 
code CLOUDY to recreate a typical narrow line region region, to derive the ionization parameter, 
and to compare our spectral data with diagnostic diagrams. 

Subject headings: active galactic nuclei — individual Mrk 78 — near infrared 



1. Introduction 

The availability of advanced near-infrared 
(NIR) instrumentation with panoramic arrays of- 
fers the possibility of employing diagnostic tech- 
niques that were formerly reserved to optical spec- 
troscopy. Besides, long-slit NIR spectroscopy of- 
fers the possibility of characterizing the velocity 
field in AGN, using different line species from high 
ionization emission lines (e.g., [Si VI] A 1.962 /im), 
H recombination lines, Fc lines (mostly triggered 
by shocks), molecular H lines, or even stellar ab- 
sorption features. These constitute a plethora of 
diagnostic tools that can be used to characterize 
many of the relevant phenomena acting in AGN, 
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from the regions close to the actual nucleus, to 
partially ionized zones, to the borders of molec- 
ular clouds, etc. Furthermore, dust extinction 
is reduced in the IR with respect to the optical 
range, thus revealing lines whose optical counter- 
parts could be obscured. 

The velocity dispersions measured in the nar- 
row emission line region (NLR) of Seyfert galaxies 
are believed to be governed by the galaxy's grav- 
itational field (Whittle 1992; Nelson & Whittle 
1996). However, in a minority of Seyferts, with 
relatively luminous linear radio sources, there is 
evidence for significant additional (jet related) ac- 
celeration (Ferruit 2002; Veilleux et al. 2002). 
The interplay between radio-emitting flows and 
the line-emitting gas is important in the study of 
active galaxies, especially those with observable 
radio jets (Cecil et al. 2000, 2002; Ferruit et al. 
1999). This interplay can occur in a wide range of 
circumstances, for instance, radio jets can be dis- 
rupted or deflected by the material in the nuclear 



1 



and interstellar media. In turn, the jets can ablate, 
shock accelerate, or destroy line-emitting clouds, 
or even trigger bursts of star formation. Mrk 78 
is a classical example of an AGN with a jet-gas 
interaction, showing large bipolar flows (Whittle 
ct al. 1988; Pedlar et al. 1989), with about half 
of the NLR flux originating in two components, 
having velocities of 680 km s -1 and —480 km s _1 , 
respectively. 

In this paper we present NIR long-slit data of 
Mrk 78 with the aim of understanding and charac- 
terizing the mechanisms responsible for its emis- 
sion spectrum. Mrk 78 is a Seyfert 2 with galac- 
tic extinction amounting to E(B — V) = 0.035 
mag (Schlegel et al. 1998). Its measured red- 
shift (z = 0.0371) (Michel & Huchra 1988) gives 
a distance of 150 Mpc and a physical scale of 715 
pc/arcsec (H = 75 km s _1 is used throughout this 
paper). HST WFPC and FOC images taken in 
optical broad-band filters and [O III] A 5007 emis- 
sion (Capetti et al. 1994; Whittle & Wilson 2004) 
show a complex structure of ionized gas aligned 
with the radio axis at PA = 84°, plus a dust lane 
running across the galaxy center at PA = 135° 
(see Fig.l). No signs of broad emission line region 
(BLR) emission has been detected from spectropo- 
larimetry (Tran 1995). In the soft X-ray images, 
taken by URl/ROSAT, Mrk 78 appears extended 
and its spectrum can be fitted by a thermal plasma 
plus galactic extinction (Levenson et al. 2001). 

Section 2 describes the observations, and Sec- 
tion 3 presents the main results, covering first the 
nuclear emission to continue with the extended 
emission spectra and the kinematics obtained from 
the line profiles and line displacements. Section 
4 discusses the extinction towards Mrk 78, and 
makes an attempt at disentangling the various ion- 
izing mechanisms at play in Mrk 78. Finally, Sec- 
tion 5 summarizes the main conclusions. 

2. Observations and Data Reduction 

NIR spectra in the range 0.8-2.4 /im were ob- 
tained on the night of 2005 March 22 at the 
4.2 m William Herschel Telescope (WHT) us- 
ing LIRIS, a recently commissioned NIR cam- 
era/spectrometer (Manchado et al. 2004; Acosta- 
Pulido et al. 2003). LIRIS is equipped with a 
Rockwell Hawaii 1024x1024 HgCdTe array detec- 
tor. The spatial scale is 0.25 arcsec/pixel and 



the slit width used during the observations was 
1", allowing a spectral resolution of 600 km s _1 
and 650 km s _1 in the ZJ and HK spectra, re- 
spectively The slit was oriented along PA = 
80°, centred on the galaxy nucleus (see Fig. 1). 
Weather conditions were relatively good, although 
with sparse cirrus. The seeing during our obser- 
vations varied bewteen 1 and 1'.'2, measured from 
the FWHM of the comparison stars in our NIR 
spectra. 
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Fig. 1. — HST(WFPC) [OIII]A 5007 image of Mrk 78, 
with the position of the LIRIS slit superimposed. The po- 
sition of the brightest IR spot is marked with a diamond. 
The origin of the coordinate system is the same as reported 
Whittle & Wilson (2004) 

Observations were performed following an 
ABBA telescope-nodding pattern, placing the 
source in two positions along the slit, separated by 
15". Individual frames were taken with integra- 
tion times of 300 s and total on-source integration 
times of 60 min in each of the ZJ and HK spectra. 
The wavelength calibration was provided by obser- 
vations of argon and xenon lamps available in the 
calibration unit at the A&G box of the telescope. 
In order to obtain the telluric correction and the 
flux calibration, the nearby HD 61232 AO V star 
was observed with the same configuration as the 
galaxy. The data were reduced following standard 
procedures for NIR spectroscopy, using IRAF 1 

IRAF is distributed by the National Optical Astronomy 
Observatories, which arc operated by the Association of 
Universities for the Research in Astronomy, Inc., under co- 
operative agreement with the National science Foundation. 
http://iraf.noao.edu/ 



2 



and the LIRIS-QL dedicated software. Consec- 
utive pairs of AB two-dimensional spectra were 
subtracted to remove the sky background. The 
resulting frames were then wavelength-calibrated 
and flat-fielded before registering and coadding all 
frames to provide the final spectrum. The final 
wavelength-calibrated galaxy spectra were divided 
by a composite spectrum to remove telluric con- 
tamination. This composite spectrum was gen- 
erated from the observed spectra of AO V star 
HD 61232, divided by a Vega model convolved 
with the actual spectral resolution as described 
by Vacca et al. (2003). Differences in the strength 
of telluric features likely due to mismatch of air 
masses and variation of atmospheric conditions 
between observations of the galaxy and the ref- 
erence star are taken into account using Beer's 
law. The IRAF task telluric was used in this step. 
The flux calibration was carried out by normal- 
izing with the JHK magnitudes provided in the 
2MASS survey catalogue (Cutri et al. 2003). The 
agreement in the continuum flux level in the over- 
lap region for the ZJ and HK spectra is quite 
good in spectral shape and absolute value (about 
15%). 

3. Results 

3.1. Nuclear Emission 

In order to study the nuclear emission we have 
extracted a spectrum covering 1.25 arcsec centred 
on the maximum of the galaxy profile. We have 
identified the maximum of the continuum as the 
position of the active nucleus and registered in the 
[O III] maps (see Fig. 1), as presented in Whittle 
& Wilson (2004). These authors have done the 
registration of the optical (with HST spatial res- 
olution) and radio images using absolute ground- 
based astrometry (Clements 1981). They found 
a good alignment between the peak of the optical 
images and the position of the radio core. This 
fact gives confidence about the NIR peak should 
also coincide with the position of their nucleus. 
The resulting nuclear spectra in the Z J, H and K 
ranges are plotted in Figs 2, 3, and 4, where the 
wavelength has been translated to the observer's 
rest frame. The most obvious features in the nu- 
clear spectra arc [S III]AA0.907, 0.953, He IA1.083, 
[Fe II]AA1.256, 1.643, Pa/3, Paa, and the coronal 
line [Si VI] A 1.962 (all wavelengths are given in 



lira). The recombination line Br7 is also present 
in the K band, although it appears very weak 
and is probably immersed in an absorption fea- 
ture. For this reason its measured flux is subject 
to a large uncertainty. It is noteworthy to mention 
that in other Seyfert 2 galaxies, e.g., NGC 2273 
and NGC 4569, this line is not found (Rhee & 
Larkin 2005). 

Emission from the H2 molecule is clearly de- 
tected, whose most prominent lines are the transi- 
tions H 2 1-0 S(2) and H 2 1-0 S(l). Both H 2 2-1 S(l) 
and H 2 1-0 S(0) emission lines are detected at < 2a 
level. 

Several other lower-intensity lines are detected 
in the ZJ spectra, namely [Ca I] A 0.985, [S VIII] A 
0.991, He IIA 1.012, [S II] A 1.032 (a blend of four [S 
II] transitions) and [P II] A 1.188, which has partic- 
ular relevance as a discriminator for different ex- 
citation mechanisms of the gas (Oliva et al. 2001). 

Wherever possible all features were measured 
by fitting a Gaussian component using the pro- 
gram twofitlines (Acosta-Pulido 2000), developed 
in the IRAF enviroment, and the resulting values 
are reported in Table 1. The [Si VI] and H 2 1 
OS (3) lines are strongly blended and two Gaus- 
sians were employed to separate the emission from 
these lines (the Starlink program DIPSO was used 
to perform the fits). The resulting fluxes and EWs 
are also reported in Table 1. The relative inten- 
sity between the H 2 1-0 S (3) and H 2 1-0S(1) lines 
is 0.83, which is in good agreement with the the- 
oretically predicted value of 1 for a wide range of 
conditions (Draine & Woods 1990), hence encour- 
aging confidence concerning the deblending proce- 
dure. The [Si VI] line appears broader (~ 70 A or 
1100 km s _1 after deconvolution) in comparison 
to other detected infrared lines (~ 600 km s _1 ). 
Rodnguez-Ardila et al. (2002) claimed a trend of 
increasing line widths with increasing ionization 
potential, perhaps indicating the existence of an 
emission region for the coronal lines closer to the 
nucleus than the classical NLR. 

3.1.1. Absorption-Line Features: The Stellar 
Population 

Several absorption features are readily seen in 
the H band nuclear spectrum of Mrk 78 (see Fig. 
3). In particular we observed the following ones: 
Mg I 1.50, CO 1.58, Si I 1.59, CO 1.62, CO 1.64 
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Fig. 2. — Flux-calibrated spectrum corresponding to the nuclear region of Mrk 78 within an aperture of 1.25", taken in the 
ZJ range. Note the simultaneous detections of high and low ionization lines, characteristic of AGN-likc spectra. 
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Fig. 3. — Same as Fig. 2 but in the H band. Note the conspicuous absorption features. This indicates a strong stellar 
contribution. 
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and Mg I 1.71 (all wavelengths are given in /im). 
These features are named according to their main 
contributor, although, depending on the stellar 
spectral type, other species may become domi- 
nant. For example, for very cool stars, OH domi- 
nates the Si I feature (Origlia et al. 1993; Dallier 
et al. 1996). In the K band, we have also detected 
Ca I 2.26, Mg I 2.28 and the 12 CO(2,0) band 
head, at 2.29 fim. Unfortunately, the entire band 
is not included in our observations as part of it 
lies outside our spectral range due to the galaxy 
redshift. 

The ratio of the equivalent width of the 1.62 
/jm CO (6,3) feature to that of the 1.59 /im Si 
I one is a good temperature indicator for late- 
type stars (Origlia et al. 1993; Forster Schrcibcr 
2000; Ivanov et al. 2004), since the CO(6,3) fea- 
ture grows rapidly from early K to late M stars, 
while the Sil one is only weakly dependent on stel- 
lar temperature. Moreover, dilution and redden- 
ing effects are cancelled out, due to the closeness 
of the features in wavelength. 

In order to identify the spectral types that pro- 
duce the absorption features seen in the HK range 



of Mrk 78, we compare our spectra with digi- 
tally available stellar templates from Dallier et al. 
(1996), observed in the H band with a medium res- 
olution. The equivalent width of the intrinsic stel- 
lar features have been computed after convolution 
of the templates, in order to match the resolution 
of our galaxy spectra. As a result of this compar- 
ison, for Mrk 78, the equivalent-width ratio of the 
CO 1.62 feature to Si 1 1.59 corresponds to spectral 
types in the range K5 III to M3 III (see Table 2). 
This is consistent with results obtained by Oliva 
et al. (1995) for a range of galaxy types and with 
those of Thatte et al. (1997) for NGC 1068. An 
age of several hundred Myr is estimated for the 
stellar system, based on the fact that intermedi- 
ate to low mass giants dominate the integrated IR 
luminosity (Renzini & Buzzoni 1986). 

In Seyfert galactic nuclei, these late-type stel- 
lar features are substantially diluted by non- 
stellar nuclear emission. Assuming a single 
stellar population, the dilution fraction (D = 
1 — EWobs/EWi n t) can be computed from the ra- 
tio of measured to intrinsic EWs of any absorption 
feature (Oliva et al. 1995). Thus in the H band, 
using the CO 1.62 feature, a starlight dilution 
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Table 1 : Observed nuclear emision lines fluxes and 
EWs. The H 2 1-0S(0) and H 2 2-1S(1) fluxes re- 
ported here correspond to an upper limit of 2a. 



Line Nucleus A Flux EW 

(pm) (lCT^erg-cm- 2 ^- 1 ) (A) 

21.5 
42.5 
3.03 
1.62 
3.26 
4.05 
31.3 
2.38 
5.93 
10.9 
6.25 
40.5 
8.34 
17.6 
3.56 
6.88 
4.13 



in the range 30-37% is obtained for the nuclear 
spectrum. For the K band CO(2,0) 2.29 we use, 
as intrinsic EWs for giant stars, those computed 
from the digitally available spectra from Wallace 
& Hinklc (1997), resulting in a dilution factor of 
43-47% in this band (see Table 2). The stellar con- 
tribution in the J band was estimated from the 
colors of K-M giants. The remaining non-stellar 
spectrum resembles a power law (F v oc v^ 1 - 5 ). 
The index of this power law is relatively close to 
the average value (F v oc i/~ 1,35 ) claimed for the 
NIR pure AGN emission (Neugebauer et al. 1979), 
which is steeper than in the optical range. 

Another common way of determining the age of 
the stellar population in the NIR is by means of 
the strength of the CO absorption feature starting 
at 2.29 /im. There are several working definitions 
of the CO index. We have used that given by 
Ivanov et al. (2000), which is more appropriate to 
the CO band measured in Mrk 78, where the spec- 
trum ends beyond 2.31 /xm. Indeed, in Mrk 78, be- 



cause of its redshift, the CO band feature is close 
to the edge of the atmospheric window; hence, 
measuring the continuum level on both sides of 
the feature is subject to large uncertainties. The 
Ivanov CO index is narrower than that of Doyon 
et al. (1994) and has the additional advantages of 
being insensitive to extinction, and to possible un- 
certainties in the continuum shape of the infrared 
spectra. We measure it as follows: 

CO = -2.51og((F 2 . 295 )/<F 2 . 282 )), 

where {F2.2&2) and (^^95) are the averaged 
flux within a bandwidth of 0.01 /im, centred on the 
blue continuum and at the band head respectively. 
The resultant CO index is 0.149 ± 0.002, which 
is among the lowest values obtained for Seyfert 
2 galaxies by Ivanov et al. (2000), as expected 
because these authors use a larger spatial scale. 
Our value is also lower than the average one cor- 
responding to the sample of pure starburst galax- 
ies as reported by the same authors. Knowing the 
dilution factor from the non-stellar component in 
the K band, which is approximately 45%, we can 
estimate the CO index after removing the non- 
stellar contribution, being its value 0.289 ± 0.007. 
In order to estimate an age from this CO index we 
need to transform its value into an age-calibrated 
one. First we convert our narrow spectral index 
into the broader, most usual one, using an expres- 
sion from Ivanov et al. (2000): 

CO(Ivanov) = (0.97±0.03)-CO(Doyon)+(0.002±0.006) 

We then obtain a photometric index using a 
linear relationship from Doyon et al. (1994): 

CO spcc - 1.46 • CO pho t - 0.02, 

which gives a CO p hot of 0.216 ± 0.023. Looking 
at Fig. 1 in Origlia & Oliva (2000), it is ambiguous 
as to which isochrone is most applicable to deter- 
mine the age of the stellar population. It is well 
known that the CO band strength is sensitive to 
surface gravity and depends both on metallicity 
and effective temperature. For solar metallicity, 
the value of the index points to ~30 Myr using 
Geneva tracks (Schaller et al. 1992), or to two dif- 
ferent values, namely ~25 and ~80 Myr, if Padua 
tracks (Bertelli et al. 1994) are chosen. These am- 
biguities have led to the use of the CO index as an 
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age indicator being questioned. Origlia & Oliva 
(2000) claim that an inadequate treatment of the 
AGB phase produces CO features that are too 
weak, so that the results are fairly dependent on 
the evolutionary tracks one adopts. Nevertheless, 
if late-type giants were the dominant population 
in the nucleus of Mrk 78, ~100 Myr would be a 
lower limit to the age of the stellar population, 
and the ~80 Myr point would then be the only 
one possible for this case. Summarizing, we con- 
clude that stellar features in the H and K bands 
are dominated by intermediate to low mass giants, 
with an estimated aged of approximately 100 Myr. 
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Fig. 5. — H band spectra of Mrk 78 (top) and two giant 
stars (K5 III and M3 III types). The spectra of the stars 
have been convolved and rebinned in order to obtain the 
same resolution as that of the galaxy spectrum. All spectra 
are normalized and shifted vertically for clarity. 



3.2. Extended Narrow Line Region 

Extended emission on both sides of the nu- 
cleus is observed, along the slit, in several spectral 
lines; namely, in [Fc II] A A 1.256,1.643, [S III] A A 
0.907,0.953, He IA 1.083, Paa, Pa/3 and [Si VI] A 
1.962. The emission in weaker lines such as Br7 
plus in several H2 molecular lines also appears 
extended although it becomes difficult to mea- 
sure, as they are at about one a level. The spec- 
tral features were measured by means of Gaussian 
fits using the in-house software utility twofitlines 
(Acosta-Pulido 2000) . The spectra were extracted 
after coadding three columns along the spatial di- 
rection, equivalent to Of! 75. However, in a region 



around 1'.'5 east of the nucleus several lines are 
split into two kinematic components (discussed in 
Section 3.3). In those cases the line fluxes are ob- 
tained by simple integration over the line profile. 

Several line flux profiles are shown in Fig. 6, 
where the Paa profile has been chosen as a base- 
line against which to compare the other line pro- 
files. We note that the flux distribution of Paa 
is asymmetric around the nucleus. It shows first 
a slow decline within l'/5 east of the nucleus and 
then falls steeply beyond 2". Instead, on the oppo- 
site side of the nucleus the profile falls quickly 1" 
west of the nucleus, further showing a shoulder at 
<~ 3". The Paa emission profile can be compared 
with the high spatial resolution Mrk 78 [O III] and 
radio images from Whittle & Wilson (2004). The 
radio image shows two lobes on both sides of the 
central nucleus, which is actually the brightest ra- 
dio source. However, in [O III] radiation there are 
two cones of emission east and west of the cen- 
ter (Whittle & Wilson 2004), which is consistent 
_ with the existence of a nuclear obscuring struc- 
ture. The large size of the [O III] obscured region 
suggests that obscuring material is required on a 
much larger scale in addition to the few parsec 

- size torus (Pier & Krolik 1992) obscuring the ac- 

- tive nucleus (Capetti et al. 1994). A transverse 
1 64 [O III] profile reveals that the emission is brighter 

l'/5 to the east, decreasing towards the west after 
2", just where the Paa shows a shoulder. 

The high-ionization line ratio [Si VI] /Paa 
shows its maximum at the position of the nucleus. 
It also shows an enhancement 2'.' 5 east of the nu- 
cleus, while it stays low (~ 0.1) in the western 
lobe. This is a clear indication that the ionization 
is higher in the nucleus and also in the east lobe, 
compared to the west side. 

In contrast, the [S III] /Paa line ratio shows a 
minimum at the nucleus and then increases to a 
nearly constant value in the lobes. Both these 
facts arc indicative of a decrease in the ionization 
level outside the nucleus, as discussed below in 
Section 4.3. 

The [Fe II] /Paa line ratio has an interesting 
behavior, showing a minimum at the position of 
the nucleus, while a high value, exceeding by a 
factor ~ 2 the value in the eastern lobe, appears 
in the western one. 

The Pa/3/Paa ratio is an indicative of extinc- 
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Table 2: Absorption line feature Equivalent Widths in the H an K bands for Mrk 78 and for three giant stars 
classes. (*M0III star) 
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tion: there are some indications that the gas in the 
east and west lobes does not suffer so much extinc- 
tion as in the nucleus, as concluded by Whittle & 
Wilson (2004). 

The He IA1.083/Paa ratio shows a shallow peak 
coincident with the nucleus of the galaxy, indicat- 
ing higher ionization than in the outer nebular re- 
gions. 

Finally, the line ratios including H 2 indicate 
that the molecular emission is compact and con- 
centrated towards the nucleus, although we cannot 
exclude the possibility that the molecular emis- 
sion may extend in a direction perpendicular to 
the ionization cone as mentioned by Reunanen et 
al. (2003). 

In summary, the scenario that we may envis- 
age through these ratios consists of a prominent 
highly ionized nucleus, which is obscured by dust; 
an eastern lobe where the ionization parameter is 
also very high according to the [Si VI] emission (a 
line with a very high ionization potential); and a 
western lobe where the [Fe II] emission is mostly 
due to the higher radio-jet interaction in the west. 
Both lobes are less affected by extinction than the 
nucleus, as can be seen from the behavior of the 
Pa/3/Paa ratio in Figure 6. This coincides with 
the existence of a dust band crossing the galaxy 
center, as can be seen in Figure 1. 

3.2.1. Eastern and Western Emission 

By looking at the variation of the line widths 
along the slit, we noted that there are two regions 
where the FWHM increases noticeably (see the 
discussion in Section 3.3). These positions coin- 
cide with the shoulders seen in the Paa flux pro- 
file. For this reason we decided to study these 
zones separately. We extracted two spectra in the 
ZJ and HK ranges, one centred on l'/5 to the 
east, and the other 2" to the west, covering 1.5 



arcsec (see Figs 7, 8, and 9). 

In the east lobe, the most prominent lines (He 
IA1.083, [Fe IIJA1.257, Pa/3 and Paa) appeared 
double-peaked, as can be seen in the insets of 
Figures 7 and 9. The two kinematic components 
correspond to a blue one, which is present overall 
along the slit, plus a weak red component, which is 
only detected close to this lobe. In the west zone, 
the FWHM appears slightly larger than the mean 
value, although it is not large enough to detect 
more than one kinematic component. 

We report integrated fluxes for the red and blue 
components of the eastern and western regions 
in Table 3. The fluxes were measured by Gaus- 
sian fitting as described in Section 3.2. In the 
case of the east region we employed two compo- 
nents, the velocity difference between them being 
~ 600 km s _1 , whereas in the west region a single 
component was used. 

The relative intensity of [Fe II] emission appears 
much brighter in the west; the values at the nu- 
cleus and at the east are similar. [Si VI] shows 
the largest value at the nucleus, and the lowest 
in the west lobe; intermediate and similar values 
are measured in the two components at the east 
lobe. Differences between the blue and red spectra 
on the east side are clearly seen in the [S III] and 
Pa/3 fluxes relative to Paa (see Table 3), showing 
a lower value for the blue component. 

The depth of the stellar features appears very 
different in the east and west lobes (see Figs 8 
and 9), being much larger in the east one. Quan- 
titatively, in the east the dilution factor is around 
20%, and in the west this factor is larger than 60%, 
measured from the CO(2,0) band head. Similar 
values are determined from the stellar absorption 
features in the H band. As expected for regions 
outside the nucleus, in the east lobe the stellar 
contribution dominates the continuum emission. 
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Fig. 6. — The most prominent emission line flux profiles vs distance from the nucleus, divided by the Pact contribution. The 
dashed line represents the flux distribution of the Pact line on a logarithmic scale. 
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In contrast, in the western lobe there could be an 
intense non-stellar continuum diluting the stellar 
counterpart, which could be related to the blue 
extended continuum claimed by Whittle & Wilson 
(2004), although its origin is still unclear. Other- 
wise, there could be an enhancement of old stellar 
population, that would contribute less to absorp- 
tion features. 




Fig. 7. — Flux-calibrated spectrum corresponding to l'/5 
to the east (top spectrum) and 2'.' to the west (bottom spec- 
trum) within an aperture of 1.5 arcsec, in the ZJ range. 
The inset shows an ampliation of [FcII]A1.257 and Pa/3 
double-peaked profiles in the east lobe. 




^.(Jlm) 



Fig. 8. — Same as Fig. 7 but in the H band. Note the 
difference in the depth of stellar features between the east 
and west lobes. 




Fig. 9. — Same as Fig. 7 but in the K band. The inset 
shows an ampliation of Pact double-peaked profile in the 
east lobe. 



3.3. Kinematics of the Ionized Gas 

A visual inspection of the two-dimensional 
spectrum reveals that the most intense emission 
lines describe what is most readily explained as a 
rotation curve, redshifted towards the east of the 
galaxy and blucshiftcd towards the west. More- 
over, all lines appear to be resolved with an aver- 
age line width of 500 km s _1 . 

In order to analyze the velocity field we have 
fitted Gaussians to the line profiles, as described 
in the previous section. We note that around 1'.'5 
in the east and near 2" in the west the measured 
line widths increase noticeably. The velocity field 
appears to be clearly perturbed in those regions. 
Furthermore, the presence of double-peaked pro- 
files becomes evident in some of the lines around 
l'/5 to the east, e.g., He IA1.083, [Fe II]A1.257, 
Pa/? and Paa (see Figs 7 and 9). We have plot- 
ted in Fig. 10 the velocity and FWHM variation 
along the slit. This has been done for an emis- 
sion line where we clearly distinguish two com- 
ponents, namely Paa, and for another line, the 
[S III] A0.953, for which a single, though broader, 
component can be fitted. In both cases, the overall 
behavior resembles that of a galaxy rotation curve 
with a half amplitude of ~ 300 km s _1 . 

The Paa line profile, in the eastern lobe, can be 
split into two kinematic components, a blueshifted 
one, roughly corresponding to the nominal rota- 
tion field, and a redshifted one, with a larger mean 
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Table 3 

Nuclear, SW and NE (red and blue components) fluxes normalized to Pack for the most 

PROMINENT LINES, ALL OF THEM FITTED BY GAUSSIANS. 



Line Nucleus 


A(/im) 


Nuclear emission 


SW emission 


Blue NE emission 


Red NE emission 


[577/] 


0.907 


0.82 ± 0.08 


1.02 ± 0.07 


0.68 ± 0.23 


1.16 ± 0.25 


[SIII] 


0.953 


1.59 ± 0.11 


2.22 ± 0.13 


1.45 ± 0.11 


2.56 ± 0.16 


Hell 


1.012 


0.14 ± 0.02 


0.12 ± 0.02 


0.13 ± 0.02 


0.13 ± 0.02 


Hcl 


1.083 


1.25 ± 0.08 


0.89 ± 0.06 


1.16 ± 0.07 


0.96 ± 0.07 


[Fell] 


1.256 


0.22 ± 0.04 


0.61 ± 0.05 


0.29 ± 0.04 


0.22 ± 0.05 


Pa/3 


1.282 


0.40 ± 0.07 


0.45 ± 0.05 


0.32 ± 0.04 


0.51 ± 0.06 


[Fell] 


1.643 


0.20 ± 0.06 


0.63 ± 0.06 


0.29 ± 0.04 


0.15 ± 0.04 


Paa 


1.875 


1.00 


1.00 


1.00 


1.00 


[SiVI] 


1.962 


0.45 ± 0.06 


0.09 ± 0.04 


0.24 ± 0.02 


0.26 ± 0.03 



velocity of ~600 km s 1 , which may correspond to 
gas accelerated by interaction with the radio jet, 
as proposed by Whittle & Wilson (2004). 

As stated above, the FWHM also increases at 
2" towards the west (see Fig. 10), suggesting that 
the gas dynamics in this region are also affected 
by the interaction with the radio jet. However, the 
spectral resolution of our spectra does not allow us 
to distinguish the different kinematic components. 



4. Discussion 

4.1. Extinction towards the Nucleus 

According to the unification scheme for Seyfert 
galaxies, the nucleus of Mrk 78 should be hidden 
from our view. This is apparently the ('MSG, clS 
for instance Capetti et al. (1994) claim strong ob- 
scuration towards the nucleus of Mrk 78, based 
on different arguments. They found that the ob- 
served UV flux is at least three orders of magni- 
tude smaller than the flux required to ionize the 
NLR. Based on this fact, applying the reddening 
curve from Cardelli et al. (1989) and the standard 
dust-to-gas ratio, they estimated a lower limit for 
the column density of Ah > 0.5 x 10 22 cm~ 2 . They 
also concluded that the nuclear emission must be 
intrinsically anisotropic in order to explain the 
observed low infrared luminosity as compared to 
the nuclear luminosity required to account for the 
overall NLR luminosity. Besides, Whittle & Wil- 
son (2004) showed color images in which the nu- 
clear extinction is clearly seen across several hun- 




Fig. 10. — Variation of velocity and FWHM (both in 
km s _1 ) measured from the lines Pact (crosses) and [S 
IIIJA0.953 (triangles). The overall measurements (fitted 
by a single Gaussian) are represented by continuous lines. 
For the Pao, which can be split into two components, 
the dashed line (diamonds) represents the blue component 
whereas the square indicates an average value for the red 
one. 
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dred parsecs. Indeed, we have not detected any 
broad wings in the H recombination lines present 
in our spectra. Previously, Goodrich et al. (1994) 
reported the non-detection of a broad component 
for Pa/3 in this galaxy, although their spectra have 
much lower signal to noise ratio than ours. 

Moreover, the extinction towards the nucleus 
can be estimated from our spectra using several 
reddening indicators, such as the [Fe II] Al .257 
and [Fe II]A1.644 lines, or other line ratios be- 
tween the observed Paschen and Brackett H re- 
combination lines. From these line ratios we 
have determined the optical extinction, Ay, us- 
ing Draine's parameterization (Draine 1989) A\ oc 
A -1 ' 75 . The results are presented in Table 4. The 
[Fe II] 1. 257/1. 644 /an line ratio is a reliable indi- 
cator given that it is not dependent on the density 
or temperature of the line-emitting gas. We mea- 
sured a value of 1.13 (sec Table 4), which yields 
a moderate value for the extinction when com- 
pared to the theoretical value of 1.34 (Bautista 
& Pradhan 1998). The [Fe II] 1. 257/1.644 line ra- 
tio measured for Mrk 78 is very close to the mean 
value 0.98 found in a sample of Seyfert galaxies 
by Rodriguez- Ardila et al. (2004). These authors 
argued that either the gas emission [Fe II] lines 
are affected by the same large extinction for most 
Seyfert galaxies, or the theoretical value is overes- 
timated, given the uncertainties in the determina- 
tion of the transition coefficients. 

Other available reddening indicators are the 
Pa/3/Paa and the Paa/Br7 line ratios. However, 
these ratios depend on the density and tempera- 
ture of the line-emitting gas. In particular, the 
Pa/3/Paa ratio varies by ~ 8 % when the tem- 
perature changes from 5000 to 20000 K (Hum- 
mer & Storey 1987). We have adopted theo- 
retical values corresponding to T e = 10 000 K 
and N e = 10 cm~ 3 . In our case, the observed 
Paa/Br7 line ratio is larger than the theoretical 
value. This could be due to an underestimation of 
the Br7 line flux, because of either the weakness 
of the line or an underlying absorption feature. 

The values derived for the visual extinction 
agree very well with those estimated by other au- 
thors using narrow emission lines (Veilleux et al. 
1997; Rodriguez- Ardila et al. 2004). It is worth 
mentioning that the extinction measured using 
these narrow lines provides only an indication of 
the material located between us and the NLR. 



However, the nucleus itself and the BLR may be 
hidden by a larger amount of material. In fact, 
the non-detection of broad wings in recombination 
lines implies a lower limit of A v > 10 (Veilleux et 
al. 1997), if we assume the existence of a hidden 
BLR. 

4.2. The Emission-Line Spectrum 

The simultaneous observation of very high and 
low ionization lines imply that a wide variety of 
physical conditions must coexist in the nuclear re- 
gion of Mrk 78. This is a characteristic feature of 
gas photoionized by a power-law radiation source 
extending from the UV through X-rays. In partic- 
ular the detection of strong coronal lines such as 
[Si VI]A1.962 (IP = 166.7 eV) and [S VIII]A0.991 
(IP = 280.9 eV) are indicative of the presence of 
extreme UV and X-ray photons (Prieto & Viegas 
2000). 

Marconi et al. (1994) concluded that the [Si 
VI] A 1.962 emission line is associated uniquely 
with Seyfert 1 and 2 nuclei, and can only be pro- 
duced by photoionization by a hard nuclear con- 
tinuum. Furthermore, they indicate that this line 
must be produced in a dense gas with large column 
densities. Hot stars, present in HII regions are not 
capable of exciting the [Si VI] line, neither increas- 
ing the ionization parameter nor the ionizing star 
temperature. This coronal line is therefore an ex- 
cellent tracer of Seyfert activity in galaxies with 
visually obscured nuclei. Indeed, Oliva & Moor- 
wood (1990) found strong [Si VI] in NGC 1068 de- 
spite the strong obscuration present in this Seyfert 
galaxy. Coronal lines are forbidden transitions of 
highly ionized species, which are formed in ex- 
treme energetic environments. Reunanen et al. 
(2003) observed that coronal line emission appears 
to be extended along the ionization cone axis, 
implying an anisotropic nuclear radiation field. 
Furthermore, coronal lines are seen with similar 
strength and frequency in both Seyfert types (Pri- 
eto & Viegas 2000), thus their study provides a 
clean test of the Unification Model. 

4-2.1. The Origin of the Intense [Fe II] Emission 

It is well known that [Fe II] emission is weak 
in HII regions but strong in the shock-excited fil- 
aments of supernova remnants. In AGN, strong 
[Fe II] emission is also common, although there 
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Table 4 

Theoretical and measured emission line ratios and A v values, calculated using Draine's 

parametrisation . 



Ratio 


Theoretical 


Measured 


Ay 


[Fell] 1.26/1.64 


1.34 


1.13±0.34 


1.97 


Pa/3/Paa 


0.49 


0.40±0.07 


1.85 


Paa/Br7 


12.2 


17.3±1.6 





is still some controversy about what is the dom- 
inant process responsible for this emission. Sev- 
eral processes may contribute to the production 
of the [Fe II] lines: 1) photoionization by soft X- 
ray-extreme UV radiation from the central source, 
producing large partially ionized regions in NLR 
clouds of high optical depth; 2) the interaction of 
radio jets with the surrounding medium, which in- 
duces shocks and produces partially ionized cool- 
ing tails; and 3) fast shocks associated with super- 
nova remnants present in starburst regions. The 
[Fe II]A1.257/Pa/3 line ratio, or equivalently [Fe 
II]A1.644/Paa, has proved to be very useful for 
distinguishing between a stellar and non-stellar 
origin for the [Fe II] emission. These line ra- 
tios increase from HII regions (photoionization) 
to supernovae remnants (shock excitation), pass- 
ing through starburst and active galaxies (Alonso- 
Herrero et al. 1997; Rodriguez- Ardila et al. 2004). 
Galaxies with [Fe II]A1.644/Paa 2 lower than 0.11 
are clasified as starbursts generally, LINERS if 
this value is larger than 0.75, and Seyfert galax- 
ies are contained in the range 0.15-0.75 (Larkin 
et al. 1998; Rodriguez- Ardila et al. 2004). In the 
case of Mrk 78, the [Fe II]A1.644/Paa line ratio 
is seen to increase outside the nucleus and shows 
larger values in the west lobe (see Section 4.3). 
We found [Fe II]A1.644/Paa = 0.2 for the nu- 
clear emission and similar values in the east lobe 
that are well within the range of Seyfert galaxies. 
Those values can be explained as due to photoion- 
ization by hard UV nuclear radiation (see Fig. 
12). Starburst activity is not sufficient to explain 
these values. In the west lobe a higher value, [Fe 

2 We transform the [Fe II] A1.257/Pa/3 line ratio to [Fe 
II] A1.644/Pact using the theoretical relations for both [Fe 
II] and H Paschen lines. The following equation was used 
[Fe II]A1.644/Paa = 0.37 x [Fe II] A1.257/Pa/3. 



II] A1.644/Paa = 0.63, was measured that is closer 
to the LINER locus, probably related to shock ex- 
citation, associated with a more efficient interac- 
tion with the radio jet. In order to explain the [Fe 
II] A1.644/Pao; line ratio in the west lobe a com- 
posite model including shock excitation plus pho- 
toionization by the nuclear radiation is required. 

On the other hand, we have detected [P II] A 
1.188 /im, which is an important emission line 
for discriminating between the different excitation 
mechanisms acting in the NLR. In particular, this 
line is very useful when compared with the [Fe II] A 
1.257 line (Oliva et al. 2001). Both lines are pro- 
duced in partially ionized regions having similar 
critical densities and excitation temperature. In 
contrast, iron is a well known refractory species 
and is strongly depleted in dust grains, whereas 
phosphorus is a non- refractory species. Photoion- 
ization alone is unable to destroy the tough iron- 
based grains, while these are easily sputtered by 
shocks. The [Fe II] /[P II] ratio is high (> 20) in 
fast shock-excited regions and low (< 2) in nor- 
mal photoionized regions (Oliva et al. 2001). For 
Mrk 78 this line ratio is 1.82 in the nuclear re- 
gion, indicating that photoionization by a soft X- 
ray continuum is the most likely mechanism re- 
sponsible for the observed nuclear spectrum. 

The Extended [Fe II] Emission and Its Relation- 
ship to Radio Emission 

As pointed out by Rodriguez- Ardila et al. (2004), 
the [Fe II]A1.644/Paa line ratio is a good discrim- 
inator to discern between a stellar or non-stellar 
origin for the [Fe II] emission, although it is not 
helpful in distinguishing between processes in the 
latter case. In addition, shocks may efficiently 
destroy dust grains and release iron into the gas 
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phase, leading to an enhancement of the [Fe II] 
emission. An excellent way to differentiate be- 
tween the different processes is the comparison of 
the [Fe II] emission and radio emission morpholo- 
gies, despite the difficulty in obtaining [Fe II] im- 
ages. Blietz et al. (1994) showed that the mor- 
phology of the [Fe II] emission correlates with that 
of the radio jet in NGC 1068. The most intense 
[Fe II] emission appears slightly offset upstream of 
the radio bow shock (see Fig. 2 in Blietz et al. 
(1994)). Mrk 78 may be a similar case in which 
the morphology of the [O III]-emitting gas seems 
to reflect different processeses in the west and east 
lobes. The east lobe can be explained as a rela- 
tively homogeneous gas structure illuminated by 
the nuclear radiation source, which is penetrated 
by the radio emission. In contrast, the west lobe 
resembles a knotty gas distribution where the in- 
teraction with the radio jet has dispersed cloud 
fragments and swept out gas, producing compres- 
sion of the gas (Whittle & Wilson 2004). 

4-. 2. 2. The Origin of the H2 Emission 

Several molecular hydrogen emission lines are 
present in the K band spectrum of Mrk 78. These 
lines are clearly detected in the nucleus and in the 
east lobe (see Figs 4 and 9), although the emission 
is mostly concentrated in the nuclear region. Rc- 
unanen et al. (2003) reported that H2 emission is 
spatially resolved preferentially in directions per- 
pendicular to the ionization cone, which is in con- 
trast with the results reported by Quillen et al. 
(1999), where H 2 emission coincides with [O III] 
and Ha emission. 

There are three possible H2 excitation mecha- 
nisms in galaxies: i) UV pumping (fluorescence), 
ii) thermal excitation due to the presence of 
shocks, and iii) hard X-ray photons, which can 
penetrate deep into molecular clouds. These three 
mechanisms produce different spectral features, 
and the relative emission line intensities can be 
used to identify the dominant mechanism. In 
particular the H 2 2-1 S(l)/l-0 S(l) line ratio is 
lower for thermal excitation (0.1-0.2) than for 
UV fluorescence (~0.55), as proposed by Mouri 
(1994). The value obtained for the nuclear re- 
gion of Mrk 78 is < 0.38, which is larger than 
the predictions for pure thermal excitation, but 
smaller than those for non-thermal UV excita- 
tion. However, in the case of dense gas (N e > 10 4 



cm -3 ), collisional de-excitation of the H2 molecule 
modifies the spectrum, approaching the thermal 
spectrum. Excitation by hard X-ray are ruled out 
by other authors based on considerations of ener- 
getics (Rodriguez-Ardila et al. 2004) and on the 
detection of the transition H 2 2-1 S(3) in several 
Seyfert galaxies (Davies et al. 2005), which would 
otherwise be suppressed in the case of X-ray irra- 
diated gas (Davies et al. 2005). 

Another way to discriminate between the ther- 
mal and fluorescent excitations is through the ro- 
tational and vibrational temperatures. In the case 
of thermal excitation both temperatures should 
be similar, whereas in the case of fluorescent ex- 
citation a high vibrational temperature will be 
in contrast with a lower rotational temperature. 
We have determined both temperatures using the 
expressions given by Reunanen et al. (2002) us- 
ing the H 2 line ratios 2-1 S(l)/l-0 S(l) and 
1-0S(0)/1-0S(2). The values obtained for this 
galaxy are T vib < 4405 K and T rot > 1557 K, 
which are in good agreement with those derived 
by Rodriguez-Ardila et al. (2004) for a number of 
Seyfert 2 galaxies, such as Mrk 279 and NGC 5728. 
Hence, the most likely mechanism for the excita- 
tion of H 2 in Mrk 78 is UV fluorescent in a dense 
gas. Note, however, that, given the large uncer- 
tainties in the line flux measurements, thermal ex- 
citation cannot be ruled out. 

4.3. Comparison with Photoionization Mod- 
els 

We have used the photoinization code CLOUDY 
(version C05.07.06), in order to reproduce the line 
ratios obtained from our spectra (for a detailed 
description see Ferland (2003)). We computed 
a grid of models based on photoionization by a 
power-law continuum and physical conditions in 
the gas typical of NLR. For the ionizing continuum 
we used a power-law shape of the form /„ ~ v a , 
where a = — 2 is the spectral index in the range 
above 10 /zm. At energies below 10 /im the spec- 
tral index is a = 5/2. We assume a plane-parallel 
geometry, metallicity equal 0.3 solar, and grains 
with properties similar to those of the Orion Neb- 
ula, which corresponds to the best-fitting models 
according to Alonso-Herrero et al. (1997). Three 
different values of hydrogen density were explored: 
n H = 10 4 , 10 5 , and 10 6 cur 3 . An input to the 
code is the ionization parameter U, defined by 
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U = Q H /(4^ 2 n H c). 

Line ratio diagrams are shown in Figs 11 and 
12 for the following combination of lines: He 
IIA1.012/Paa versus [S III]A0.953/Paa and [Si 

VI] A1.963/Paa versus [Fe II] A1.644/Paa. Dashed 
lines represent the computed line ratios for the 
three hydrogen densities employed, beginning with 
an ionization parameter of log U = and decreas- 
ing it with a step of 0.5, towards the bottom of 
these plots. 

Values of these line ratios from Table 3, cor- 
responding to the nucleus, the west and the east 
lobes (red and blue components for the latter), are 
overplotted for comparison with our simulations 
made with CLOUDY. It can be seen that data 
from our spectra are compatible with photoioniza- 
tion models for all regions considered within val- 
ues of log U in the range —1 to —2, except for 
the [Fe II]/Paa line ratio. By looking at the dia- 
gram in Fig. 11, the best fit for the nucleus and 
the east lobe-blue component is obtained when 
logU ~ —1.5. However, for the west lobe and east 
lobe-red component, the best fit is obtained when 
logU <~ —1.8. In all cases the H density is close 
to riH = 10 5 cm~ 3 , which approaches the critical 
density for the [S III] line. From Fig. 12, it is 
easy to see that the [Si VI]/Paa line ratio is very 
sensitive to the value of the ionization parameter, 
becoming saturated at values logU ~ —0.5. At 
higher values of the ionization parameter, the [Si 

VII] A2. 480 transition will become more important. 
The nuclear emission can thus be explained with 
values log U ~ —0.5, whereas in the east lobe this 
value decreases to logU <~ —1.3, and in the west 
lobe to log/7 ~ —1.6. The differences between the 
best-fitting values of U from the different line ra- 
tios may be reconciled by lowering the spectral 
index of the power law, i.e., hardening the ion- 
izing continuum, or using combination of clouds 
with different optical depths. In order to explain 
the [Fe II] /Pace ratio, models may have to include 
effects due to interaction with the radio emission, 
such as increasing the iron abundance after grain 
sputtering. 

We therefore conclude that the measured line 
ratios from the spectra of Mrk 78 are compatible 
with photoinization by an AGN-like continuum. 
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Fig. 11. — HeIIA1.012/Paa versus [SIII] A0.953/Paa di- 
agnostic diagram computed with the photoionization code 
CLOUDY. Models are calculated for a power— law contin- 
uum (index -2), different values of H density (different line 
types) and a sequence of ionization parameters. The marks 
along the U— sequence are separated by logU=0.5. Our mea- 
surements are also overplotted: cross for the nucleus; as- 
terisc for SW lobe; square for the blue component in the 
NE lobe; diamond for the red component in the NE lobe. 
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Fig. 12. — Same as Fig. 11 but with [Fell] A1.644/Paa 
versus [SiVI] A1.963/Paa The meaning of symbols is the 
same as in Fig. 11 

5. Conclusions 

We have presented and analyzed the NIR line 
spectrum of the nucleus and extended NLR of the 
Seyfert 2 galaxy Mrk 78. The following results 
were found: 

• The nuclear spectrum and the extended 
NLR spectrum are produced by photoion- 
ization by a hard UV-X-ray continuum, as 
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indicated by the presence of the [Si VI] and 
[P II] lines. This fact is in good agreement 
with the presence of a hidden active nucleus 
in Mrk 78. 

• The nuclear emission in the H band has nu- 
merous stellar absorption features. Accord- 
ing to the relative equivalent widths, the 
stellar population is dominated by K-M gi- 
ants in this wavelength range. The nuclear 
stellar emission is diluted by a non-stellar 
component by a factor ~ 34% in the H band 
and ~ 45% in the K band. 

• An enhancement of the [Fe II] emission is 
observed outside the nucleus, in particular at 
~ 2'.' to the west. This effect is interpreted 
in terms of the interaction with the radio 
emission. 

• Several transitions of the H 2 molecule are 
detected in the K band. The dominant ex- 
citation mechanism could not be determined 
unambiguously, although we favor UV fluo- 
rescence in a dense gas. 
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